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Abstract. We provide an integral representation for (implied) copulas 
in terms of the moment generating function of dependent random vari- 
ables. The proof uses ideas from Fourier methods for option pricing and 
can be applied to a large class of models from mathematical finance, 
including Levy and affine processes. 



1. Introduction 

Copulas provide a complete characterization of the dependence between 
random variables, and link in a very elegant way the joint distribution with 
the marginal distributions via Sklar's Theorem. They are however a rather 
static concept and do not blend well with stochastic processes which can be 
used to describe the random evolution of dependent quantities, e.g. the evo- 
lution of several stock prices. Therefore other methods to create dependence 
in stochastic models have been developed. Multivariate stochastic processses 
spring immediately to mind, for example Levy or affine processes (cf. e.g. 
ISato 19991 IDuffie et al. 20031 ICuchiero et al. 2011|) . while in mathematical 
finance models using time-changes or linear mixture models have been devel- 



oped; see e.g. 


Luciano and Schoutens (2006) Luciano and Semeraro (2010) 


Kawai (2009) 


Eberlein and Madan (2010) or 


Khanna and Madan (2009) 



to mention just a small part of the existing literature. In these approaches 
however the copula is typically not known explicitly. Another very interest- 



ing approach due to Kallsen and Tankov (2006) introduced Levy copulas to 



characterize the dependence of Levy processes. 

The aim of this short note is to provide a new representation for the 
(implied) copula of a multidimensional random variable in terms of its mo- 
ment generating function. The derivation of the main result borrows ideas 
from Fourier methods for option pricing, and the motivation stems from the 
knowledge of the moment generating function in most of the aforementioned 
models. This paper is organized as follows: in section [2] we provide the rep- 
resenation of the copula in terms of the moment generating function; the 
results are proved for random variables for simplicitly, while stochastic pro- 
cesses are considered as a corollary. In section [3] we provide two examples to 
showcase how this method can be applied, for example, in performing sen- 
sitivity analysis of the copula with respect to the parameters of the model. 
Finally, section [3] concludes with some remarks. 
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2. Copulas via Fourier transform methods 

Let M n denote the n-dimensional Euclidean space and (-, •) the Euclidean 
scalar product. The scalar product is extended from real to complex numbers 
as follows: for u = (iti)i<i< n and v = (wj)i<i<„ in C n , set (u, v) := Ya=i u i v i'i 
therefore we do not use the Hermitian inner product Y27=i u ^%- Moreover, 
w := (i«i)i<j< n - 

Let us consider a random variable X = (X±, . . . ,X n ) T G W 1 defined on 
a probability space (Q,, P). We denote by F the cumulative distribution 
function (cdf) of X and by / its probability density function (pdf). Let C 
denote the copula of X and c its copula density function. Analogously, let 
Fi and fi denote the cdf and pdf respectively of the marginal X{, for all 
i G {1, . . . , n}. In addition, we denote by F^ 1 the generalized inverse of Fi, 
i.e. F" 1 ^) = ini{v £ R : F^v) > u}. 

We denote by Mx the (extended) moment generating function of X: 

M x (u) = E[e^], (2.1) 

for all u G C n such that Mx(u) < oo. In the sequel, we will assume that the 
following conditions are in force. 

(HI): Assume that M X (R) < oo, for R € 1 C R n . 

(H2): Assume that M X (R + i-) G L 1 (R n ) for R G J C E n . 

(H3): Assume that K := 1 n J n \{0} ^ 0. 

Theorem 2.1. Let X be a random variable that satisfies assumptions (HI)- 
(H3). The copula of X is provided by 

I r ~{R+iv,x) 

C{u) = —— M x {R + w) dv (2.2) 

(-2vr) n J iU=i{Ri + m) x i =F i 1 (ui) 

where u G [0, l] n and R G TZ. 

Proof. Assumption (H2) implies that F\, . . . , F n are continuous and we know 
from Sklar's Theorem that 

C(u u ...,u n ) = F(F 1 -\u 1 ),...,F- 1 (u n )); (2.3) 



sec 



McNeil, Frey, and Embrechts (2005, Theorem 5.3) for a proof in this set 



ting and Riischendorf (2009) for an elegant proof in the general case. 

We will evaluate the joint cdf F using the methodology of Fourier methods 
for option pricing. That is, we will think of the cdf as the "price" of a digital 
option on several fictitious assets. Let us define the function 

9(y) = l{ 1n <x u ..., yn <x n }(v), x,y€ M n , (2.4) 
and denote by g its Fourier transform. Then we have that 
F{x)=TP(X 1 <x 1 ,...,X n <x n ) 

= ^{ 1 {X 1 <x l ,...,X n <x n }] 

= j^yi J M X {R + iv)g(\R - v)dv, (2.5) 
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where we have applied Theorem 3.2 in Eberlein, Glau, and Papapantoleon 



(|2010p . The prerequisites of this theorem are satisfied due to (H1)-(H3) and 
because gn € L 1 (R n ), where <7r(x) := e ( - R ' x ^g(x) for R e 1Z. 

Finally, the statement follows from (|2.3p and ()2.5[) once we have computed 
the Fourier transform of g. We have for Ri < 0, i S {1, . . . , n}, 



g(iR - v) 



j(iR-v,y) 



I 



{yi<xi,...,y„ 



n 

i=l 



^(-Ri-m)yi 



(-irn 



-{Ri+ivi)xi 



L (Ri + ivi) 



which concludes the proof. 



(2.6) 



□ 



Remark 2.2. If the moment generating function of the marginals is known, 
the inverse function can be easily computed numerically. We have that 

F' 1 ^) = \nf{v eR:Fi(v)> u} 

= inf{ueR:E[l {x .< w} ] >u}, 



where the expectation can be computed using (|2.5p again, while a root find- 
ing algorithm provides the infimum (using the continuity of Fi). 

We can also compute the copula density function using Fourier methods, 
which resembles the computation of Greeks in option pricing. 



Lemma 2.3. The copula density function c of X is provided by 
1 



c(u) 



J M x (R + iv)e- {R+iv ' x) dv 



Xi=F i 1 (u i )' 



(2.7) 



(2*)" IlLi /<(**) 
where u € (0, l) n and R € 7i. 

Proof. Let u € (0, l) n , then we have that X{ = F^ l (ui) is finite for every 
i € {l,...,n}, hence e~^ R ' x ^ is bounded. Using assumption (H2) we get 
that the function Mx(R + iv)e~( R+w ' x ^ is integrable and we can interchange 
differentiation and integration. Then we have that 



c(u) 



duj . . . du n 

qu 



C(ui, ...,u n ) 

1 



du\ . . . du n (— 27r) T 



M x (R + iv) 



q—(R+iv,x) 



M x (R + iv) 



(-2n) n J \\ n i=1 {R i + w i )du 1 ...du n 



-(R+iv,x) 



dv 

i =F- 1 (u l ) 

dv. (2.8) 



Xi=F i 



1 



ANTONIS PAPAPANTOLEON 



Now, using the chain rule and the inverse function theorem we get 

9 n f e - {R+iVtX) \ 



dui . . . du n \ \x i =F~ 1 ( : u l ) 

n 

f = \ Ui=ifi[ x i) *i= F i K) 

which combined with (|2.8p yields the required result. □ 

A natural application of these representations is for the calculation of cop- 
ulas of random variables from multidimensional stochastic processes. There 
are many examples of stochastic processes where the corresponding charac- 
teristic functions are known explicitly; prominent examples are Levy pro- 
cesses, self-similar additive (Sato) processes and affine processes. 

Corollary 2.4. Let X = (Xt)t>o be an M n -valued stochastic process on a 
basis ($7, (J r j)t>o,P). Assume that the random variable Xt, t > 0, satisfies 
assumptions (H1)-(H3). Then the copula of X t is provided by 

1 f e -(R+iv,x) 

°t( u ) = / M Xt {R + w) dv (2.10) 

where u 6 [0, l] n and R € 1Z. An analogous statement holds for the copula 
density function ct of X t . 

3. Examples 

We will demonstrate the applicability and flexibility of Fourier methods 
for the computation of copulas with two examples. First we consider a 2D 
normal random variable and next a 2D normal inverse Gaussian (NIG) Levy 
process. Altough the copula of the normal random variable is the well-known 
Gaussian copula, little was known about the copula of the NIG until recently; 



see Theorem 5.13 in Schmidt (2003) for a special case. Now lv. Hammersteinl 
([2011} Ch. 2) has provided a general characterization of the (implied) copula 
of the multidimensional NIG law using properties of normal mean- variance 
mixtures. 

Example 3.1. The first example is simply a "sanity check" for the proposed 
method. We consider the 2-dimensional Gaussian distribution and compute 
the corresponding copula for correlation values equal to p = {—1,0, 1}; see 
Figure 13.11 for the resulting contour plots. Of course, the copula of this 
example is the Gaussian copula, which for correlation coefficients equal to 
{— 1, 0, 1} corresponds to the countermonotonicity copula, the independence 
copula and the comonotonicity copula respectively. This is also evident from 
Figure EU 

Example 3.2. Let X = (Xt)t>o be a 2-dimensional NIG Levy process, i.e. 

X t = (X},Xt) ~ NIG 2 (a, (3, St, fit, A), t > 0. (3.1) 

The parameters satisfy: a, 5 € R+, /3, fx G M 2 , and A € M 2x2 is a symmetric, 
positive definite matrix (w.l.o.g. we can assume det(A) = 1). Moreover, 
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Figure 3.1. Contour plots of copulas for Example 13.11 



o? > (/3,A/3). The moment generating function of X±, for u € M 2 with 
a 2 - ((3 + u, A((3 + u)) > 0, is 

M Xl (u) = exp ({u, ^) + 5 (y« 2 - (J, A/3) - y/a 2 - {(3 + u, A((3 + ) , 



cf. Barndorff-Nielsen (1998) The marginals are also NIG distributed and we 
have that X\ ~ NIGtd',^,^,^), where 



a 



Si, 



, /3 J = /3 l + /3 i ,5 2 ^ 1 , * = $ 



for i = {1,2} and j = {2, 1}; cf. e.g. Blassild (1981, Theorem 1) Conditions 
(HI) and (H2) are satisfied for fi£K z such that c? - {P + R, A{/3 + R)) > 0, 



see Appendix B in Eberlein et al. (2010) , thus 1Z ^ 0. 

Therefore, we can apply Theorem 12.11 to compute the copula of the NIG 
distribution. The parameters used in the numerical example are similar to 
Eberlein et al. (2010, pp. 233-234) | a = 10.20, (3 = 



-3.80 \ § 



0.150 and 



~ V -2.50 )i 

fi = 0, and two matrices A + = ( q ? ) and A - = ( _i "2 ) , which lead to pos- 
itive and negative correlation. The correlation coefficients are /?+ = 0.1015 
and p- = —0.687 respectively. 

The contour plots are exhibited in Figures 13.21 and 13.31 and show clearly 
the influence of the different mixing matrices A + and A~ to the dependence 
structure. A minor influence of time in the dependence structure can be also 
observed, in particular on the examples with positive correlation. 



4. Final remarks 

We will not elaborate on the speed of Fourier methods compared with 
Monte Carlo methods in the multidimensional case; the interested reader is 



refered to Hurd and Zhou (2010) for a careful analysis. Moreover, Villigcr 
(|2007p provides recommendations on the efficient implementation of Fourier 
integrals using sparse grids in order to deal with the "curse of dimensional- 
ity". However, let us point out that the computation of the copula function 
will be much quicker than the computation of the copula density, since the 
integrand in (|2.2p decays much faster than the one in (|2.7p . One should think 
of the analogy to option prices and option Greeks again. Finally, it seems 
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2D NIG, positive correlation 2D NIG, negative correlation 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 '0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Figure 3.2. Contour plots of copulas for NIG (t = 1). 



2D NIG, positive correlation 2D NIG, negative correlation 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 '0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Figure 3.3. Contour plots of copulas for NIG (t = 7). 

tempting to use these formulas for the computation of tail dependence co- 
efficients. However, due to numerical instabilities at the limits, they did not 
yield any meaningful results. 
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